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ABSTRACT 



Context. X-ray spectra of many accreting pulsars exhibit significant variations as a function of flux and thus of mass accretion rate. 
In some of these pulsars, the centroid energy of the cyclotron line(s), which characterizes the magnetic field strength at the site of the 
X-ray emission, has been found to vary systematically with flux. 

Aims. GX 304-1 is a recently established cyclotron line source with a line energy around 50 keV. Since 2009, the pulsar shows 
regular outbursts with the peak flux exceeding one Crab. We analyze the INTEGRAL observations of the source during its outburst in 
January-February 2012. 

Methods. The observations covered almost the entire outburst, allowing us to measure the source's broad-band X-ray spectrum at 
different flux levels. We report on the variations in the spectral parameters with luminosity and focus on the variations in the cyclotron 
line. 

Results. The centroid energy of the line is found to be positively correlated with the luminosity. We interpret this result as a manifes- 
tation of the local sub-Eddington (sub-critical) accretion regime operating in the source. 

Key words. X-ray binaries - neutron stars - accretion 



1. Introduction 



In accreting binary pulsars, matter from the normal stellar com- 
panion is transferred to a highly magnetized (B > 10 12 G) neutron 
star. In the vicinity of the accretor, the gas flow is disrupted by 
the neutron star's magnetic field and channeled towards the mag- 
netic poles, where most of the X-rays originate. The physics and 
the structure of the X-ray emitting region(s) ab ove the neutron 
star surface are still highly debated (see, e.g., iBecker & Wolff! 
I2007t iFarinelli et alJl20ia . Since the matter hitting the accre- 
tor's surface is highly ionized, the magnetic field strength is a 
crucial parameter determining the physical processes inside the 
emitting region and the formation of the observed X-ray spec- 
trum. A direct way to assess the B-field strength at the site of X- 
ray emission is the measurement of the cyclotron resonant scat- 
tering features (CRSF or cyclotron lines) in the X-ray spectrum 
of a pulsar. These features appear as absorption lines, caused by 
the resonant sc attering of photons off the electrons in La ndau 
levels (see, e.g.. iTriimper et al1ll978t Ifsenberg et al. 1998). The 
energy of the fundamental line and the spacing between the har- 
monics are directly proportional to the field strength. 



In some accreting pulsars, the energy of the cyclotron line 
has been found to vary with luminosity, apparently due to 
a displacement of the line formation region. Such variations 
of the line ener gy hav e been reported for V 0332+53 (e.g., 
[Tsygankov eTaD 120101). 4U 01 15+63 (e.g.. Ffsygankoy et al.1 
20 10). Her X-1 (IStaubertet al.ll2007l: iKlochkov et al.ll2011h . and 
A 0535+26 (Klochkov et al. 2011). The luminosity-dependence 
of the cyclotron feature has strong implications f or the physics 
of the X-ray emitting region a s disc ussed, e.g., bylStaubert et al] 
(120071) . IKlochkov et all (l20Tll) . and lBecker & Wolffl (l2007l) ~ 

GX 304-1 is a rec ently established cyclotron line source 
dYamamoto et al.l [201 lh . It was discovered in a balloon exper- 
iment in 1967 and subse quently identified as an X-ray pulsar 
with a period of ~272s (McClin tock et alJ 1977b. The system 
contains a Be-type optical comp anion dMason et alJ [l978) and is 
located at a distance of ~2.4kpc dParkes et al . 1980). Since 1980, 
GX 304-1 has remained in a quiescent state, showing no out- 
bursts. Startin g from 2008, when the source was detected with 
INTEGRAL dManousakis et al J2008I) . GX 304-1 "resumed" its 
activity exhibiting outbursts with a period of ~ 132.5 d. 

The energy of the cyclotron line in GX 304 -1 was measured 
with Suzaku and RXTE to be around 52 keV bv lYamamoto et all 
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Table 1. INTEGRAL observations of the GX 304-1 outburst in 
January-February 2012. 



2011 



Jan 20 



Jan 31 



Feb 10 



Rev. 


Obs. ID 


Mid. MJD 


Exposure [ksec] 








JEM-X 


IBIS 


SPI 


1131 


09400230006 


55944.0 


64.6 


42.7 


68.6 


1132 


09400230007 


55947.0 


42.4 


31.9 


36.6 


1133 


09400230008 


55950.0 






10.7 


1134 


09400230009 


55952.8 


7.3 


25.4 


37.8 


1135 


09400230010 


55955.7 




6.7 


25.1 


1136 


09400230011 


55958.7 


36.9 


28.1 


32.9 


1137 


09400230012 


55962.0 


78.1 


59.7 


78.4 


1138 


09400230013 


55965.0 


60.7 


45.2 


52.3 



(1201 lh . These authors demonstrate that the data taken at different 
flux levels show an indication of a positive correlation between 
the line energy and the X-ray flux, although at a low signifi- 
cance level. This made GX 304- 1 a good target for a luminosity- 
dependent study of the cyclotron feature. In this work, we 
present the analysis of INTEGRAL data (see next section for 
the mission description) taken during the January 2012 outburst 
of the source. The data reveal a positive correlation between the 
cyclotron line energy and the source flux, as well as the variation 
in other spectral parameters during the outburst. We discuss our 
finding in the context of a model assuming that different accre- 
tion regimes can operate in a particular pulsar depending on its 
X-ray luminosity. 

2. Observations and data reduction 

At the beginni ng of January 2012, GX 3 04-1 entered an outburst 
as reported by Yamamot o et al.1 (120121) . The outburst was moni- 
tored by the Int ernational Gamma-R ay Astrophysics Laboratory 
- INTEGRAL dWinkler et al J 12001 . starting at MJD -55943.5, 
when the source flux in the ~20-80keV range was ~250mCrab, 
through the maximum of the outburst, when the source flux ex- 
ceeded one Crab, to MJD ~55965.5, when the flux dropped 
to ~100 mCrab. INTEGRAL performed one observation every 
satellite orbit (about three days), with a typical exposure of a 
few tens of kiloseconds each. In total, eight observations were 
performed. The INTEGRAL scientific payload contains three 
X-ray instrum ents: (i) the imager IBIS sensitive from ~20keV 
to a few MeV dUbertini et alj|2003t ILIsbrun et al.ll2003l) : (ii) the 
spectr ometer SPI sensitive in roughly the same energy range as 
IBIS dVedrenne et al.ll2003l) : and (iii) the X-ray monito r JEM-X 
operating between ~3 and ~35keV dLund et al.ll2003l) . TableQ] 
summarizes the INTEGRAL observations of GX 304- 1 . The in- 
creased solar activity during the observations led to the reduction 
in the exposure time and availability of the instruments as can be 
seen from the exposure columns of TableQ] 

The INTEGRAL observations are indicated in FigQ] 
which shows the Swift/BAT light curve of GX304-lfl The 
INTEGRAL monitoring has an excellent coverage of the entire 
outburst, providing a rare opportunity to follow the evolution of 
the outburst from the early rising phase to the late decay phase. 

For our analysis, we used data from the ISGRI detector of 
IBIS, which is sensitive in the 20-300 keV energy range, JEM-X, 
and SPI. Standard data processing was performed with version 9 
of the Offline Science Analysis (OSA) software. We performed 
an additional gain correction of the ISGRI energy scale based 




55940 55945 



55950 55955 
Time. MJD 



55960 55965 



1 We used the Swift/BAT transient monitor results provided by the 
Swift/BAT team 



Fig. 1. INTEGRAL observations of GX 304- 1 (horizontal bars) 
superposed on the Swift/BAT light curve of the source during its 
outburst in January-February 2012 (~0.22 units of the vertical 
axis corresponds to one Crab). 



on the background spectral lines of Tungsten (the modified re- 
sponse files based on the nearest Crab observations were used). 



3. Spectral analysis 

Figure [TJ shows that the X-ray flux of GX 304-1 during the ob- 
served part of the outburst changed by an order of magnitude. 
This allowed a detailed study of the luminosity dependence of 
the source's broad-band X-ray spectrum. X-ray pulsations with 
a period of ~274.9 s were detected in all INTEGRAL observa- 
tions. This value is roughly consistent with the known pulse pe- 
riod. A detailed timing analysis is not part of the present work 
and will be presented elsewhere. 

In all observations, the X-ray continuum could be closely 
modeled by a standard power-law/cutoff function (flux oc 
E r exp [E/Ef \,i], where E is the photon energy, F is the pho- 
ton index, and Ef id is the exponential roll-off parameter) modi- 
fied by photo-electric absorption at low energies. In addition, the 
spectra showed a cyclotron resonant scattering feature around 
~50keV in absorption, which was modeled with a multiplica- 
tive absorption line with a Gaussian optical depth profile I{E) = 
I cont (E) ■ e- G{E \ where G(E) = -r cyc /( V2lrcr cyc ) • exp[-0.5(£ - 
£cyc) 2 /o"c yc ], ^cont(£) is the continuum function, E cyc , cr cyc , and 
r cyc are the centroid energy, width, and optical depth of the line, 
respectively. The line is clearly detected in ISGRI and SPI data 
separately, as shown by the residual plots in Fig. [2] The inclusion 
of the absorption line in the model leads to an improvement in 
the reduced^ 2 from 3.10 for 145 d.o.f. (with the corresponding 
null-hypothesis probability of only ~10~ 32 ) to 0.75 for 142 d.o.f. 
(null-hypothesis probability >0.9). The energy of the line is con- 
sistent with that reported by I Yamamoto et al.l (1201 ll) based on 
the RXTE and Suzaku observations. We also included an addi- 
tive Gaussian component to model the Fe K a fluorescence emis- 
sion line around 6.4 keV. The inclusion of the line reduces the 
residuals around 6 ke V. The corresponding improvement of the 
reduced^ 2 is, however, marginal: from, 0.80 (143 d.o.f.) to 0.75 
(142 d.o.f.). The presence of the Fe line is, therefore, question- 
able. 

Both the continuum and cyclotron line parameters vary sys- 
tematically during the outburst. Here, we focus on the evolution 
of the cyclotron line energy E cyc , to establish a possible correla- 
tion o f E cyc with flux, similar to that found by Yamamot o et al.l 
d201 lh . To characterize this variability, we used only observa- 
tions where data from all three INTEGRAL X-ray instruments 
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Fig. 2. Example of the INTEGRAL spectrum of GX 304- 1 (rev- 
olution 1132) with a simultaneous fit of data from JEM-X, 
ISGRI, and SPI with a power-law/cutoff model including a cy- 
clotron absorption line (see text) (a), the residuals for a fit with 
the model without the cyclotron line for JEM-X + ISGRI (b), 
and JEM-X + SPI (c), and with the model where the line is in- 
cluded for data from all instruments (d). 
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Fig. 3. Evolution of the cyclotron line centroid energy E cyc 
throughout the outburst of GX 304-1 as measured with 
INTEGRAL. The vertical error bars indicate lcr-uncertainties. 
The horizontal error bars indicate the time intervals of the obser- 
vations. The dotted line represents the re-scaled Swift/BAT light 
curve. 



were available. This filtering allowed us to obtain a homoge- 
neous set of data and, hence minimize possible systematic ef- 
fects. Following this selection, data from the orbits 1133 and 
1 135 were excluded. The JEM-X data collected during the orbit 
1134 (peak of the outburst) have been flagged as "bad" by the 
instrument team owing to the impossibility to provide a precise 
energy calibration. In our analysis, however, we use the JEM- 
X data to restrict the low-energy part of the broad-band X-ray 
continuum of the source, for which a precise energy scale is not 
very critical. Therefore, after consulting the JEM-X team, we re- 
introduced the JEM-X data from the revolution 1 134 in our anal- 
ysis. We checked, however, that inclusion of the JEM-X data of 
rev. 1 134 to the corresponding spectral fit does not significantly 
influence the measured energy of the cyclotron line (the focus of 
this work), but only affects its uncertainty. Table[2] summarizes 
the best-fit spectral parameters achieved in each observations. 




-9.0 -8.5 -8.0 

logio(F 4 _ 80 kev), erg cm -2 s -1 

Fig. 4. Cyclotron line centroid energy E cyc as a function of the 
logarithm of flux in the 4-80 keV range. The error bars indicate 
lcr-uncertainties (the flux uncertainties are smaller than the sym- 
bol size). The top X-axis represents the corresponding isotropic 
source luminosity assuming a distance of 2.4 kpc. The dotted 
line shows a linear fit to the £ , cyc -logio(Flux) dependence. 



A clear systematic variation in the line energy over the out- 
burst is evident in Fig. [3] which shows that E cyc generally fol- 
lows the X-ray flux. To assess the interdependence of the two 
parameters, we plotted E cyc as a function of the X-ray flux in the 
4-80 ke V range measured with INTEGRAL in the respective ob- 
servations (Fig.|4]l. The plot shows a positive correlation between 
the two values. A linear fit to the dependence of E cyc on the loga- 
rithm of flux (dotted line in the plot) reveals a slope of 4.97+1.12 
keV/logio(erg crrT 2 s _1 ). The standard linear correlation analysis 
of the £' cyc -logio(Flux) dependence yields a Pearson's correla- 
tion coefficient of 0.88 with a probability of obtaining the corre- 
lation by chance of ~0.01 (one-sided). 

We verified whether our spectral fits contain artificial (model- 
driven) dependences of the cyclotron line energy E cyc and other 
model parameters using ^-contour plots. No significant model- 
driven dependences between E cyc and any of the continuum pa- 
rameters were found. The line energy was, however, found to be 
somewhat coupled to the line width <x cyc and its central optical 
depth r cyc . Nevertheless, the ^-minima and the confidence inter- 
vals could be clearly identified and are separated for the different 
observations as shown in Fig. [5] The plot shows the contours for 
the E cy J<r cyc pair. The E cy jT cyc contours look similar. 

To check whether the E cy Jfixix correlation is related to in- 
strumental effects, we performed spectral fits using only the 
JEM-X and SPI data (excluding ISGRI). To verify whether the 
correlation depends on the choice of the spectral model, we fit 
the data using a Lorentzian line profile instead of a Gaussian 
one. We also tried two alternative continuum functions: XSPEC 
powerlawxhighecut and compTT models. In the formermodel, 
the highecut component controls the exponential roll-off. In ad- 
dition to Efoid, this component has an additional parameter - the 
cutoff energy E cuto s, above which the spectrum is affected by the 
roll-off. In our fits, E cu toff stays between a few and ~10keV. In 
all cases, the positive E cyc /flux correlation was reproduced. We 
conclude, therefore, that the reported correlation arises from the 
source's behavior and reflects real physics. 



4. Discussion and conclusions 

The positive correlation between the cyclotron line centroid en- 
ergy and the flux found with INTEGRAL confirms the claim 
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Table 2. Best-fit spectral parameters with the corresponding 1 cr-uncertainties for the INTEGRAL observations used in this work. 
The luminosity is provided assuming a distance of 2.4 kpc and an isotropic emission diagram. 



Revolution 


1131 


1132 


1134 


1136 


1137 


1138 


r 

-Efoid, keV 
cr cyc , keV 

Flux/10 -8 erg s -1 cirr 2 
L x /10 37 erg s- 1 


1 10 +0M 

1,lu -0.03 

19 r 07 

52 98+ 80 
6 45 +0 - 62 

"■^-0.58 

9.4 + _\i 

n ^71+0.002 
v.j i J_ 002 

0.26 


0.83!°$ 
15.5!^ 
54 39+ - 75 

7 56+ - 59 
12 1 +1 ' 4 

n 8 1 1 +0.004 
U - 8i 1 -0.004 

0.56 


n qo+o.09 

"• 7J -0.09 

55.28tS 
8 45 +a83 
13 3 +24 

1 630 +00U 

i.uju _ oou 
1.12 


16.91™ 
51 44+0.54 

1 -0.52 

5 60 +0 - 48 

•'•""-0.46 

e+0.9 
8.8_ 09 

O 79Q+0.005 
"' ' ^-0.005 

0.50 


i 90 +0.08 

i - i ->-0.08 
19 5 +L1 

51 70+ - 66 

'"-0.63 
« 44+0.56 
J ' -0.55 

8 9 +L2 
-1.1 

340+ 002 

u.jtu _ 003 
0.23 


1 56 +0 13 

i - JU -0.05 

28.6 +3 ; 8 
48.37 + |-]« 
4.801™ 

i5r 0002 
j -0.002 

0.10 


xL/d.o.f. 


1.1/164 


0.7/164 


1.3/164 


1.1/164 


0.9/164 


1.0/164 



> 8- 



b 6 - 




Fig. 5. x -contour plots of the parameter pair E cyc /cr C y C for a 
few selected observations. The contours correspond to ^ in + 1.0 
(the projections of this contour to the parameter axes correspond 
to the 68%-uncertainty for one parameter of interest), ^ in +2.3 
(68%-uncertainty for two parameters of interest), and ;t^ in +4.61 
(90%-uncertainty for two parameters of interest). The respective 
orbit numbers are indicated. 



a radiative shock. As discussed in Staub ert et alj ([2007) and 
iBecker et al.l d2012l) . the height of the emitting region decreases 
with increasing luminosity owing to a corresponding increase 
in ram pressure of the infalling material, leading to a posi- 
tive Ecyc/flux correlation. The critical luminosity L c depends 
on the parameters of the accreting n eutron star, but is gener- 
ally around a few times ~10 37 ergs~ 1 dBasko & Sun vaevll 19761 
Stau bert et al.l2007l:lBecker et al.ll2012h . Assuming a distance of 
2.4 kpc dParkes et al.lll980l) . the X-rav luminosity of GX 304-1 
in the 4-80 keV range during the reported INTEGRAL observa- 
tions varies between ~l.lxl0 36 erg s _1 and ~1.13xl0 37 ergs _1 . 
Thus, according to the described picture, the source should be- 
long to the class of "sub-critical" sources, for which a positive 
Ecyc/flux correlation is expected. The reported observations are, 
therefore, in agreement with the idea of two accretion regimes 
and increases the yet very small sample of accreting pulsars with 
established E cyc /flux correlations. 
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of this correlation by Yamamot o et alj (1201 lb . The revealed 
dependence establishes GX 304-1 as the third member of a 
slowly emerging class of accreting pulsars showing a positive 
E^vc/flux correlation with the other members being HerX-1 
(Stau bert et all 120071) and possibly A 0535+26 dKlochkov et alJ 
1201 ll) . for which the positive correlation has so far been only 
established in the pulse-to-pulse analysis. The opposite (neg- 
ative) correlation between the line energy a nd the flux was 
found in V 0332+53 and 4U 01 15+63 (e.g.. iTsvgankovet alJ 
20071120 lOl). Accordin g to discussions inlStaubert et al.l (l2007l) ~ 



Klochkov et alJ d201 lh . and lBecker et alJ d2Q 1 2h . these two types 
of dependences reflect two different regimes of accretion. A 
particular regime is realized in a source depending on whether 
its X-ray luminosity L is above or below a critical luminos- 
ity L c , which corresponds to the local Eddington luminosity at 
the X-ray emitting structure(s) on/above the neutron star sur- 
face. In accreting pulsars radiating above L c ("super-critical" 
sources), infalling matter is decelerated in a radiative shock, 
whose height is believed to increase with L, i.e., drift towards 
an area with a lower B-field strength. The opposite behav- 
ior probably occurs in sources radiating below or close to L c 
("sub-critical" sources), where infalling matter is stopped by 
the Coulomb drag and collective plasma effects rather than in 
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